Enhancing the spreading of quantum walks on star graphs by additional bonds 
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We study the dynamics of continuous-time quantum walks (CTQW) on networks with highly 
degenerate eigenvalue spectra of the corresponding connectivity matrices. In particular, we consider 
the two cases of a star graph and of a complete graph, both having one highly degenerate eigenvalue, 
while displaying different topologies. While the CTQW spreading over the network - in terms of 
the average probability to return or to stay at an initially excited node - is in both cases very 
slow, also when compared to the corresponding classical continuous-time random walk (CTRW), 
we show how the spreading is enhanced by randomly adding bonds to the star graph or removing 
bonds from the complete graph. Then, the spreading of the excitations may become very fast, 
even outperforming the corresponding CTRW. Our numerical results suggest that the maximal 
spreading is reached halfway between the star graph and the complete graph. We further show how 
this disorder-enhanced spreading is related to the networks' eigenvalues. 

PACS numbers: 05.60.Gg, 05.60.Cd, 
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I. INTRODUCTION 

Transfer processes play a very important role in many 
physical, chemical or biological instances, involving en- 
ergy, mass or charge. Such transport processes strongly 
depend on the topology of the system under study, as, for 
instance in simple crystals [21], complex molecules [10] or 
spin networks [3]. 

One often encounters regular networks, such as lat- 
tices or rings, where most or all of the nodes have the 
same degree. Examples are the star graphs (SG), and the 
complete graphs (CG), which have attracted lately much 
attention in the field of quantum computation. Such 
studies focused on the quantum central limit theorem for 
continuous-time quantum walks on SG [19], the discrete- 
time Grover walk on a SG with one loop [11], and quan- 
tum searches on highly symmetric (complete) graphs [18]. 
Exact analytical results were recently obtained for quan- 
tum walks on SG and on CG whose Hamiltonians have 
one highly degenerate eigenvalue. This, as also discussed 
in what follows, may let the quantum walk remain local- 
ized at its starting node, which then renders the quantum 
spreading slow. 

Spreading of an excitation can be made faster by ran- 
domly adding (deleting) bonds to the SG (from the 
CG). Here we do not allow the additional bonds to con- 
nect one node with itself and forbid double bonds be- 
tween any pair of nodes. The inclusion of the additional 
bonds changes the structure of the network and, conse- 
quently, the eigenvalue set of the corresponding Hamil- 
tonian. Furthermore, destroying the original regularity 
of the SG (CG), the additional bonds create "shortcuts" 
in the graph, by which the spreading through CTQW 
may occur faster, as, for instance, found for small-world 
networks [16]. 

For modelling transport processes on networks one 



can use several approaches: in quantum mechanics, the 
Hamiltonian is determined by the connectivity of the sys- 
tem. For example, the dynamics of an electron in a crys- 
tal can be described by the Bloch ansatz [21]. This can 
be also related to transport processes in polymers, where 
the connectivity plays an important role in the dynamics 
[6]. 

Classical transport processes modelled by continuous- 
time random walks (CTRW) are described by a master 
equation, involving a transfer operator T based on the 
topology of the system [10], [8], [9]. A quantum me- 
chanical analogue of CTRW, namely, one variant of the 
continuous-time quantum walk (CTQW), can be intro- 
duced by using the transfer operator T in defining the 
Hamiltonian H [7]. For simple lattices this is equiva- 
lent to a nearest- neighbor hopping model [7], [12], [13], 
[17]. The transformation replaces the classical diffusion 
process by a quantal propagation through the structure. 

The paper is organized as follows. Section II introduces 
the CTRW and CTQW concepts. Tools for determining 
their spreading are displayed in Section HI. Section IV fo- 
cuses on CTRW and CTQW on SG and on CG, whereas 
the spreading over graphs intermediate between the SG 
and the CG is studied in Section V. Section VI concludes 
with a summary of results. 



II. QUANTUM WALKS ON NETWORKS 



In the following we focus on the dynamics of excita- 
tions over networks; these consist of N nodes connected 
by bonds. The information about a network's topology, 
i.e., its connectivity, is stored in the N x N connectivity 
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matrix A, whose elements A^j are: 

Ajcj = ^ — 1 if fc and j are directly connected 
otherwise. 



(1) 



Here fj is the number of bonds emanating from j, also 
referred to as the degree of j. The connectivity matrix 
has the following properties: 

1. A is symmetric and real. 

2. All its eigenvalues En are real and non- negative. 

If the network is simply connected, A has a single van- 
ishing eigenvalue i?min = 0. In order to model the dy- 
namics in the two extreme cases of a purely coherent 
and of a purely incoherent (diffusive) transport, we em- 
ploy the concepts of CTQW and of CTRW, respectively. 
In both cases, the dynamics is largely influenced by the 
network's topology, i.e., by A. Now, we depict an ex- 
citation localized at node j through the state the 
ensemble of the \ j) states forms an orthonormal basis set 
— l,-.-,-^'^}- In the incoherent CTRW case, the 
spreading occurs based on the transition rates jkj j which 
denote the probability to go per unit time from j to k. 
Assuming these rates to be the same for all the nodes, 
i.e., taking j^j = 7 (and we will set 7=1 without loss 
of generality), induces a simple relation between T and 
A, namely, T = -A [8], [7], [15]. 

In the quantum case, the states \j) span the whole 
accessible Hilbert space. The time evolution of an exci- 
tation initially placed at node |j) is determined by the 
system's Hamiltonian H; in the approach of Ref. [7] 
used here, the CTQW is directly related to T through 
H = — T [7], [15]. The classical and quantum mechani- 
cal probabilities to be in state at time t when starting 
at t = from state |j) arc then: 



and 



Pfej(i) = (fc|exp(Tt)|j) (2) 
^fc,,(t) = |(fc|exp(-im)| (3) 



respectively, where we set in Eq. (3) h = 1. For t = 
both expressions read Pk,j{(^) = T^k,j{0) = Skj, where Skj 
is Dirac's delta-function. 

To get from the formal Eqs. (2) and (3) to the explicit 
solutions for a particular lattice one has to diagonalize 
T and H. Given that here the CTQW Hamiltonian H 
is the negative of the CTRW transfer matrix T, lets the 
eigenvalues and eigenvectors of both operators be practi- 
cally identical. The difference in the dynamics of CTQW 
and of CTRW is due to the different functional forms of 
Eq. (2) and of Eq. (3). 



III. CLASSICAL AND QUANTUM SPREADING 

While one needs to calculate the eigenvectors (which 
can become tedious for large networks) for making defi- 
nite statements on the transition probabilities pfej (t) and 



TTkjit), there are quantities - to be defined below - which 
only depend on the eigenvalues, whose distribution, the 
density of states (DOS) or spectral density, is given by: 



1 ^ 



(4) 



The DOS contains information about the system and 
allows to analyse several features which depend on the 
network's topology. For CTRW one such feature is the 
classical average probability to be (return or remain) at 
the initially excited node averaged over all possible initial 
nodes J [1], [2]: 



j n=l 



(5) 



Hence p{t) depends only on the eigenvalues £"„ of A. but 
not on its eigenvectors |0„). 

For CTQW, on the other hand, n{t) = ^ Y.^^i ^j.jit) 
depends on the eigenvectors \4>n)- However, one can de- 
fine, in a fashion similar to the classical p(t) , the quantity 

[12], [14], 



1 ^ 1 " 



-iE„t 



(6) 



Here a{t) also depends only on the eigenvalues En- 
Furthermore, |Q!(i)P obeys the Cauchy-Schwarz incqual- 
ity [12]: 



{t) 



> 



(7) 



Hence |a(t)p is a lower bound to 7f(t). We will use 
Eqs. (5) and (6) to assess the spreading [14]. Now, clas- 
sically, a quick decrease of p{t) means a quick increase in 
the probability of finding the excitation away from the 
initially excited node. We hence infer that the spreading 
is faster when p{t) decreases more quickly. In the quan- 
tum case, because of the unitary time evolution, 7f(t) 
and |a(t)P show oscillations. However, the overall (aver- 
age) decay of p{t) and of |Q:(i)P can be used to infer the 
spreading, say, by focusing on the decay of the envelope 
of, in particular, |Q:(t)P [14]. 

For regular d— dimensional networks these procedures 
to assess the spreading have been used and discussed in 
[14]. The result was that quantum walks over such net- 
works appear to be faster than the classical ones: the 
envelope of |Q;(t)|^ turns out to decay as i"'', whereas 
p{t) decays as t~ . Evidently, other networks may (and 
do) behave differently. In such cases we will use the long- 



3 



time averages 



RW 



1 

lim ^ j dt Pit) 

T— i-oo 1 J 



and 



P< 



QW 



lim — 

T->oo T 



dt \a{t)\ 



(8) 



(9) 



in order to extract a time-averaged global spreading mea- 
sure. 

For CTRW, due to the eigenvalue Ei = E„,in = 0, p{t) 
will eventually drop to the equipartition value as 
one can see by rewriting Eq. (5) as: 



Pit) 



1 

'n^ N 



1 ^ 

ri=2 



(10) 



It hence follows that \imT-^oo p{t) = i/N and also that 
-Prw = 1/-^^- Thus, in the incoherent case at long times 
both quantities do not depend on the topology of the con- 
sidered network. In the quantum case, |a(t)p oscillates. 
With Eqs. (6) and (9) it follows: 



QW 



= lim - /dt V 

■n.n' 



-i{E„-E„,)t 



N ^ 



p{En) 



(11) 



In Eq. (11) 5{En — En') is unity if En = En' and vanishes 
otherwise. Note that Y^Er, /'(^n) = 1- Thus, the DOS 
completely determines the long-time average of \a(t)\'^. 
For instance, for a ring with an odd number of nodes all 
eigenvalues are non-degenerate, i. e. there are N differ- 
ent eigenvalues, each with p{En) = 1/A^, and therefore 
in that case ^q'^ = 1/iV. The DOS determines hence 
the CTQW spreading measure Pqw- 
We can rewrite Eqs. (5) and (6) as 



P{t) 



E 



p{En)e 



-E„t 



and 



nit) > 



E 

E„ 



p{En)e' 



-iEnt 



(12) 



(13) 



respectively. In general, for a single highly degenerate 
eigenvalue £'m, with p{E,n) of order 0(1), while the other 
eigenvalues have a DOS at most of order 0{1/N), the 
average transition amplitude is [14]: 



iit) = [p{En.)e-^^-' + PiEn)e-'^"'}, (14) 

E„^E„ 



from which one gets to order 0(1/ N"^) the approximate 
expression: 



i{t)\^ ^ p{E„,){p{E,n)+ E P{En)cOs[{E,n-En)t]y 
E„^E„ 



(15) 

Equation (15) shows that for highly degenerate eigenval- 
ues the lower bound |Q:(t)P wiU not decay to zero, but 
will oscillate about a finite value. This means that for 
CTQW there is a high probability to remain or to return 
to the initially excited node. This leads to a slow aver- 
age spreading for CTQW on such networks. However, for 
lattices of high dimensions, there is a significant proba- 
bility that a walker never returns to the origin, thus, the 
CTQW can get transient [5]. In the following we will 
focus on two particular examples of such networks and 
on how to overcome the slow spreading. 



IV. STAR GRAPH AND COMPLETE GRAPH 

There exist many networks with highly degenerate 
eigenvalues. Here, we focus on two specific networks 
which have both one highly degenerate eigenvalue but 
vastly different topologies: On the one hand we consider 
the SG, which is a network consisting of N nodes with a 
central node connected to iV — 1 leaf nodes, which are not 
connected to each other (see Fig. 1). Therefore, the cen- 
tral node has degree A^ — 1, and each of the leaf nodes has 
degree 1. On the other hand we consider the CO, where 
each of the N nodes is connected to all other — 1 nodes 
(see Fig. 1). Both networks are very "ordered", in the 
sense that there is an "exchange" symmetry, meaning 
that exchanging the positions of any pair of nodes (ex- 
cept for the central node of the SG) leaves the network 
invariant. 




FIG. 1: SG and CG of size iV = 10 

The connectivity matrices of the two graphs read 

... _i\ 



AsG = 



f N -1 -1 ■■ 
-1 1 

■■ 







(16) 



\ -1 •■•0 1 / 
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and 



^CG 



( N -I -1 

-1 iV- 1 

: -1 



(17) 



\ -1 -1 •■• -1 A^- 1 y 

Consequently, the Hamiltonians can be written as 



N 



HsG = {N- 1)|1)(1| + ^ (b-)01 - 11)01 - |J>(1|) (18) 



and 



H 



CG^{N~l)Y,{\3){j\) 



(19) 



The eigenvahics of both graphs can be calculated ana- 
lytically. The SG of size N has only three distinct eigen- 
values [14]: Ei^Q, E2 = ■■■ = En-1 = 1, and En ^ N 
and their degeneracies are D(Q) = 1, D{1) ~ N — 2, 
and D{N) = 1, respectively. The CG of size N has 
only two distinct eigenvalues [20], namely Ei = and 
E2 = ■ ■ ■ = Em = N , with degeneracies D(0) = 1 and 
D{N) =N -I. 

The dynamics of CTQW and CTRW on SG and on CG 
has been studied in [20], where exact analytical results 
for the transition and the return probabilities have been 
obtained. It has been shown that when the central node 
in the SG is initially excited, the spreading turns out to 
be equivalent to that over a CG of the same size. 

Using Eq. (12) and (13) one obtains the following ana- 
lytical expressions for p{t) and for |Q:(t)|^ [15] for the SG 

Pit) = 1 { 1 + (iV - 2)e-* + e-(^-2)* } (20) 

and 

|a(t)|2 = l^\l + {N- 2)e-'* + e'^^^-^)* I (21) 

For the CG the analytical expressions for p{t) and 
\a{t)\^ read [15]: 



PW = ^{l + (A^-l)e-^*} 



and 



\a{t)\ 



1 



1 + (iV - l)e 



-iNt 



(22) 



(23) 



It is a simple matter to calculate the long-time aver- 
ages of the expressions above. For Eqs. (20) and (22) 
the CTRW long-time averages are Prw = 1/-^ both for 
the SG and for the CG. In the quantum case we use 
D{En) — Np{En) in Eq. (11). Such that, using now the 



explicit eigenvalues and their degeneracies listed above, 
we obtain: 



and 



pSG _ 
-^QW ~ 




-4N - 


h6 


> 


1 




iV2 




N 


pCG _ 




- 2N - 


h2 


> 


1 




iV2 




N 



(24) 
(25) 



We conclude that a quantum walker on such networks 
in particular for large N, has a large probability to stay or 
to return to the initially excited node, much larger than 
for a classical walker. We may even see this as a sign 
that the quantum spreading is slower than its classical 
analogue. Clearly, this is due to the fact that one eigen- 
value of H is in both cases highly degenerate. One may 
compare Eqs. (24) and (25) to the situation of a quantum 
walk over a ring with an odd number of nodes, for which 
^QW ~ 1/-^- In the latter case there are no degener- 
ate eigenvalues, so that D{En) = 1 for all n. Note that 
from the value -Pq^ — 1/N one should not infer that the 
CTQW is never faster than the CTRW. The result only 
implies that in the long-time limit the return probability 
is in average equally distributed among all nodes of the 
network. 

We now turn to the question in how far the situation 
changes when we (randomly) add bonds to the SG. By 
this we will explore the transition from the SG to the CG. 
As it will turn out, the CTQW spreading may increase 
when the disorder gets larger. 



V. FROM THE STAR GRAPH TO THE 
COMPLETE GRAPH 

In this section we consider graphs generated from a SG 
of size N to which B additional bonds are added. The 
new bonds are not allowed to connect any node with 
itself; furthermore only a single bond between any pair 
of nodes is permitted. The total number of additional 
bonds, Bniiix{N), needed to convert the SG into the CG 
of the same size is: 



(N - 1){N - 2) 



(26) 



Examples of SG with B = 3 additional bonds are given 
in Figure 2. 

As seen from Fig. 2, randomly adding bonds to a 
SG may lead to distinct topologies. While for B very 
close to 1 or to iV the number of topologically distinct 
graphs is small, this number increases rapidly when leav- 
ing these regions (B/N and (B - N)/N <€. I). Now, 
each topological realization leads to a set of eigenvalues. 
Graphs with distinct topologies may have an identical set 
of eigenvalues (the graphs are then called co-spectral [4]). 
For computational reasons we will not distinguish be- 
tween such graphs and wc will denote by the term "con- 
figuration" the set of graphs leading to the same set of 
eigenvalues. 





FIG. 2: Examples of SG with _B = 3 additional bonds 



We will call the set of all distinct configurations ob- 
tained from a SG by adding B additional bonds a "clan" , 
and denote the number of elements in the clan by Afs- 
Evidently, Mb depends on the size of the network. We 
determine for every fixed B value the number of distinct 
eigenvalue sets inside the corresponding clan by consid- 
ering all the corresponding Hamiltonians. This is done 
as follows: 

1. Generate the Hamiltonian H of the given SG (see 
Eqs. (16) and (18)). There arc {N - 1){N - 2)/2 
zero entries above (below) the diagonal (first row 
and column have only non-zero entries). 

2. Consider all the possibilities of changing pairwise 
such entries from Hy- — Hji = to = Hji = 
— 1, corresponding to the insertion of an addi- 
tional bond between the nodes i and j. There are 



(iV-l)(jV-2) 



ways of doing this, leading to the same 



number of distinct H. 

3. For each H obtained in this way determine the 
eigenvalue set and sort the eigenvalues in ascending 
order. 

4. Group together all the H corresponding to a fixed 
B value and determine by comparison of the as- 
cending eigenvalues the number of distinct eigen- 
value sets Afs ■ 

For the numerical procedures we use the standard 
package of Mathematica 8.0. To fix the ideas we start 
exemplarily with graphs of iV = 10 and with the corre- 
sponding clans. We arc interested in how many config- 
urations A/b arc included in the clan B and how many 
graphs are represented by a given configuration. For this 
we show in Tabic I for several B values the corresponding 
Mb- 

Moreover, it turns out that a configuration containing 
graphs which are less symmetric with respect to the ex- 
change symmetry is more probable to show up when dis- 
tributing the additional B bonds randomly. To exemplify 



this we also show in Table I the most probable and the 
least probable graph topology together with their eigen- 
value sets. One remarkable feature of Table I is that the 
least probable configurations have eigenvalue sets con- 
sisting of integer values only, while the eigenvalue sets 
corresponding to the most probable configurations also 
include non-integer polynomial roots. 

Figure 3 shows Mb as a function of B for different N. 
As can be seen from the figure, Mb reaches a maximum 
around B{N) = i?max(^)/2. One may note the very 
large increase of Mb with growing N (note the logarith- 
mic scale on the y-axis). For TV = 10 (where Bmax = 36) 
and B taken to be i3niax/2 = 18 we find more than 
3- 10'* distinct eigenvalue sets, while for N = 7 (where 
Bmax = 15) and B = 8 there are only 215. For exceed- 
ing 10 it gets quickly impossible to determine through 
the above algorithm all the configurations. Thus, for 
A > 10 we will rely on randomly creating networks and 
obtain the ensemble average using a relatively large num- 
ber R = 10000 of realizations. Formally we set: 



R 

R 



(27) 



where r runs over the particular realizations. In this way 
we determine the ensemble-averaged probabilities (p(t)) r 
and (|a(i)p)ii', along with the long-time average {Pqw)r- 
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FIG. 3: Number of distinct eigenvalue sets as a function 
of B for the networks of A^ = 7, A = 8, A^ = 9, and 
A^ = 10 nodes. Note the logarithmic scale. 



A. Average probability of being at the origin 

The infiuence of additional bonds on the spreading is 
captured in the behavior of the transition probabilities 
between certain pairs of nodes. However, as can be an- 
ticipated, the functional form of these transition proba- 
bilities strongly depends on the particular nodes chosen. 
Therefore, we start by considering the ensemble averages 
of the average probability to be (return or stay) at the 
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B 




The least 
probable 
configuration 


The most 
probable 
configuration 


The least 
probable 
eigenvalue 
set 


The most 
probable 
eigenvalue 
set 


4 


11 








{10, 5, 3, 3, 1, 1, 1, 1, 1, 0} 


{10, 4.414, 3, 3, 1.586,1, 1, 1, 1, 0} 


10 


1470 






{10, 6, 6, 6, 6, 1, 1, 1, 1, 0} 


{10, 6.26, 5.4, 4.601, 3.4743, 2.625, 2.5858, 
1.705, 1.34, 0} 


18 


31566 


^^^^ 




{10, 7, 7, 7, 7, 7, 4, 4, 1, 0} 


{10, 8.33, 7.5, 6.395, 5.6087, 5.3913, 4.605, 
3.4562, 2.66, 0} 


26 


1470 








{10, 10, 10, 10, 10, 5, 5, 5, 5, 0} 


{10, 10, 9.408, 8.8192, 8.353, 7, 6.24, 5.396, 
4.7835, 0} 










32 


11 






{10, 10, 10, 10, 10, 10, 8, 8, 6, 0} 


{10, 10, 10, 10, 10, 9.4142, 8, 8, 6.5858, 0} 



TABLE I: For networks of size = 10 and different B we display the least and the most probable eigenvalue sets. 



initial node. In particular, we focus on the ensemble av- 
erage of for CTRW and of the lower bound |Q;(t)P for 
CTQW, see Eqs. (5) and (6). Taking the time-average of 
(|Q;(i)p)fl yields {Pqw)r, see Eq. (9). Here and in the 
following we focus on graphs with = 10. 

Figure 4 shows (Pqw) r. a-s a function of the number of 
additional bonds. Clearly, one observes a decrease from 
the SG- value to a broad minimal plateau centered around 
B = i3niax/2 (see the inset), after which the values in- 
crease again until the CG-valuc is reached. The mini- 
mal value of (Pqvv)_r at B = -Bmax/2 approaches 1/A^ 
from above. Therefore, when the number of additional 
bonds is so large as to yield the maximal number of possi- 
ble configurations, the time-averaged spreading becomes 
comparable to that for the corresponding GTRW. From 
this we can already infer that adding (removing) bonds to 
(from) the SG (CG) enhances the spreading of CTQW. 
This behavior differs from the one found for CTQW on 
SWN which are build up from a ring-like configuration; 
for SWN the spreading becomes slow with increasing B 
[16]. 

In order to highlight better the behavior of the cor- 
responding CTRW, we have to take a closer look at the 
temporal development of the average probability to be at 
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FIG. 4: {Pqw)f/. for networks of A' = 10 nodes with 
R = 10000. 



the initial node. For the SG and the CG there exist exact 
analytic expressions for |a(i)P and for p{t), see Eqs. (20) 
and (23). In both cases the lower bound |a(t)P shows 
oscillations but no decay to a value comparable to 1/A^ 
(the long-time value for p{t)). Specifically, |a(i)P oscil- 
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latcs around the values given by Eqs. (24) and (25) for 
the SG and for the CG, respectively. 

Now, on adding bonds to the SG (|a(t)|^)_R follows a 
complex course: First it decreases to lower values un- 
til reaching B = i3max/2 = 18; then it increases again 
towards the behavior found for the CG. This is exem- 
plified in Fig. 5 which shows (|Q:(t)p)/j for networks of 
size TV = 10 with 5 = 4, 18, and 32 additional bonds. 
One notes that {\(jt{t)\'^) ji for B = A explores roughly 
the same interval as does (p(t))fl, see Fig. 5(a). Taking 
B = Bniax/2 = 18 yields a {\a{t)\'^)R whose time-average 
is close to Prw = 1/N but which also - for short times 
- reaches values which are about one order of magnitude 
below those of the CTRW, see Fig. 5(b). Adding more 
bonds reverses the behavior, i.e., (Pqw)^ becomes larger 
and the amplitude of the oscillations around this value 
becomes smaller, see Fig. 5(c). 
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FIG. 5: Averaged probabilities {p{t))ji and {TT{t))ii on 
double logarithmic scales for the SG of size = 10 
with B = 4 (a), 18 (b) and 32 (c), where R = 10000. 

Since {\a{t)\'^)]i is a lower bound, care has to be taken 
when interpreting the above results. It still might be 
possible that the exact value, namely {wit)) a, lies above 
the CTRW curve. However, as has been shown earlier 
[14], [15], the maxima of {\a{t)\'^)n tend to be of similar 
value as the maxima of (7f(t)}fl. In particular, for B = 
^max/2 there is a considerable short-time part (up to t w 



4) where {\a{t)\'^)ji lies below {p{t))ii and which includes 
several local maxima. Note that at t ~ 4 {p{t))fi has 
already reached the equipartition value. We thus infer 
that in this case during the (short) time interval needed 
by the CTRW to reach equipartition, the spreading of 
CTQW is faster than that of the corresponding CTRW. 

Eigenvalue sets 

Both {\a{t)\^)R for CTQW as well as (p{t))R for 
CTRW depend only on the corresponding eigenvalues (of 
the Hamiltonian and of the transfer matrix, respectively). 
As mentioned above, the SG and the CG have eigenvalue 
sets with one highly degenerate eigenvalue. This situa- 
tion will change as a function of B. 

In order to visualize the transition from the SG to the 
CG in the domain of the eigenvalues, we consider as a 
function of B the quantity: 

dB{E) = 0{E - Ei^y))R ^]{Ze{E^ [Ei^\). 

n n,r 

(28) 

Here 9{x) is the Heaviside function and En are the 
eigenvalues of the rth realization of the distribution of 
the B additional bonds. The quantity dsiE) gives the 
average number of eigenvalues below E. Figure 6 shows 
dB{E) for 5 = 4, 10, 18, 26, and 32. 
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FIG. 6: Aveage number of eigenvalues for A = 10 below 
E] here the number of additional bonds is 
B — 4:, 10, 18, 26, and 32, see text for details. 

With increasing B the eigenvalues move to higher val- 
ues, as is evident from the fact that the dsiE) curves 
decrease with increasing B. Sharp steps in dsiE) indi- 
cate eigenvalues with large degeneracy; such instances are 
particularly clear for B = A and for B = 32 whereas the 
curve for i? = 18 is much smoother. This is in line with 
the behavior of |a(t)p: a few highly degenerate eigenval- 
ues lead to a plateau in |Q!(t)P whereas non-degenerate 
eigenvalues let |a(t)p drop to values close to the CTRW 
equipartition value 1/A^. 
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VI. CONCLUSION 

We have studied the dynamics of continuous-time 
quantum walks on networks whose connectivity matrices 
have one highly degenerate eigenvalue. This fact leads to 
a slow spreading of CTQW in term of the average prob- 
ability to return or remain at the initially excited node 
of the network. In particular, we focused on two types 
of networks, the star graph and the complete graph. We 
studied the crossover from the star graph to the complete 
graph by randomly adding bonds to the star graph un- 
til reaching the complete graph. In so doing we found 
that in the ensemble average the spreading gets faster 
when adding approximately B = B^^y^jl bonds to the 
star graph. The reason for this behavior is to be found 



in the ensemble averaged distribution of the eigenvalues 
of the connectivity matrices. Adding (extracting) more 
and more bonds to (from) the star (complete) graph re- 
sults in broader distributions. Thus, we have shown that 
under disorder, obtained by cither adding bonds to the 
star graph or by deleting bonds from the complete graph, 
the quantum walks spread more. 
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